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Abstract 

This is the second in a series of three papers in which we study a two-dimensional lattice gas 
consisting of two types of particles subject to Kawasaki dynamics at low temperature in a large 
finite box with an open boundary. Each pair of particles occupying neighboring sites has a negative 
binding energy provided their types are different, while each particle has a positive activation energy 
that depends on its type. There is no binding energy between particles of the same type. At the 
boundary of the box particles are created and annihilated in a way that represents the presence 
of an infinite gas reservoir. We start the dynamics from the empty box and are interested in 
the transition time to the full box. This transition is triggered by a critical droplet appearing 
somewhere in the box. 

In the first paper we identified the parameter range for which the system is metastable, showed 
that the first entrance distribution on the set of critical droplets is uniform, computed the expected 
transition time up to and including a multiplicative factor of order one, and proved that the 
nucleation time divided by its expectation is exponentially distributed, all in the limit of low 
temperature. These results were proved under three hypotheses, and involve three model- dependent 
quantities: the energy, the shape and the number of critical droplets. In the second paper we prove 
the first and the second hypothesis and identify the energy of critical droplets. In the third paper 
we settle the rest. 

Both the second and the third paper deal with understanding the geometric properties of 
subcritical, critical and supercritical droplets, which are crucial in determining the metastable 
behavior of the system, as explained in the first paper. The geometry turns out to be considerably 
more complex than for Kawasaki dynamics with one type of particle, for which an extensive 
literature exists. The main motivation behind our work is to understand metastability of multi- 
type particle systems. 
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1 Introduction 



Section [TTT] defines the model, Section rOl introduces basic notation. Section [1751 states the main theo- 
rems, while Section [1.41 discusses the main theorems and provides further perspectives. 

1.1 Lattice gas subject to Kawasaki dynamics 

Let A C be a large box centered at the origin (later it will be convenient to choose A rhombus- 
shaped). Let 

d-A^{xeA: 3yiK: \y-x\^l}, 
9+A = {a;^A: 3yeA: |y-x| = l}, 

be the internal, respectively, external boundary of A, and put A^ = A\c'^A and A+ = A U 9+A. 
With each site a; S A we associate a variable r]{x) G {0,1,2} indicating the absence of a particle or 
the presence of a particle of type 1 or type 2. A configuration rj = {ri{x) : cc S A} is an element of 
X — {0, 1, 2}^^. To each configuration 77 we associate an energy given by the Hamiltonian 

(x,j/)eA*.- xeK xeK 

where A*'^ — {(x, y): x,y E A^, |x — j/| = 1} is the set of non-oriented bonds inside A^, —U < is 
the binding energy between neighboring particles of different types inside A^, and Ai > and A2 > 
are the activation energies of particles of type 1, respectively, 2 inside A. W.l.o.g. we will assume that 

Ai<A2. (1.3) 

The Gibbs measure associated with H is 

fip{v)^^e~f'"^^\ rye A-, (1.4) 

where /3 e (0, 00) is the inverse temperature and Zp is the normalizing partition sum. 

Kawasaki dynamics is the continuous-time Markov process, {rit)t>a with state space X whose tran- 
sition rates are 

cp{r^,T^') ^ e-f'^'''^^"^-"^^'^^^ ,7,r;'eA',,;^77',r,0 77', (1.5) 
where rj ^ rj' means that ?/ can be obtained from rj by one of the following moves: 

• interchanging and 1 or and 2 between two neighboring sites in A 
("hopping of particles in A"), 

• changing to 1 or to 2 in 9^ A 

("creation of particles in 9~A"), 

• changing 1 to or 2 to in d^A 

("annihilation of particles in d~A"), 

and cp{r], rj') = otherwise. Note that this dynamics preserves particles in A, but allows particles to be 
created and annihilated in d~A. Think of the latter as describing particles entering and exiting A along 
non-oriented bonds between 9+A and d~A (the rates of these moves are associated with the bonds 
rather than with the sites). The pairs (?/,?]') with rj -H- rf are called communicating configurations^ the 
transitions between them are called allowed moves. Note that particles in d~ A do not interact: the 
interaction only works in A^. 

The dynamics defined by (jl.2p and (jl.Sp models the behavior inside A of a lattice gas in Z^, 
consisting of two types of particles subject to random hopping with hard-core repulsion and with 
binding between different neighboring types. We may think of Z^\A as an infinite reservoir that keeps 
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the particle densities fixed at pi — e'" and p2 = e^" In the above model this reservoir is replaced 
by an open boundary d~ A, where particles are created and annihilated at a rate that matches these 
densities. Thus, the dynamics is a finite-state Markov process, ergodic and reversible with respect to 
the Gibbs measure /i^ in (|1.4p . 

Note that there is no binding energy between neighboring particles of the same type. Consequently, 
the model does not reduce to Kawasaki dynamics for one type of particle when Ai = A2. 

1.2 Notation 

To state our main theorems in Section 11.31 we need some notation. 
Definition 1.1 (a) D is the configuration where A is empty. 

(b) W is the set consisting of the two configurations where A is filled with the largest possible checkerboard 
droplet such that all particles of type 2 are surrounded by particles of type 1 . 

(c) Lo: rj ^ rj' is any path of allowed moves from r] € X to r]' G X . 

(d) ${r],r]') is the communication height between ri,r]' € X defined by 

<l>{ri,ri') ^ min maxi7(^), (1.6) 

and B) is its extension to non-empty sets A, B d X defined by 

<P(A,B)^ min ^{rj.r]'). (1.7) 

(e) Vri is the stability level of rj ^ X defined by 

V^=-P{V,Ir,)-H{Tj), (1.8) 

where T^i = ^ X : i?(^) < H(ri)} is the set of configurations with energy lower than 77. 

(f) Xstah = {rj € X : H{ri) = min^^x H{S^)} is the set of stable configurations, i.e., the set of configu- 
rations with mininal energy. 

(g) '^mota = {'76 X: Vrj = m.a,yi^izx\x^i^i,y{\ fn^toLstable configurations, i.e., the set of 
non-stable configurations with maximal stability level. 

(h) r = Vj-f for rj G A'nicta (note that r] 1-^ Vq is constant on X^cta.), T* — ^(□,ffl) — H{D) (note that 
HiD) - 0;. 

In [3j we were interested in the transition of the Kawasaki dynamics from □ to ffl in the limit as 
/? — >■ 00. This transition, which is viewed as a crossover from a "gas phase" to a "liquid phase", is 
triggered by the appearance of a critical droplet somewhere in A. The critical droplets form a subset of 
the set of configurations realizing the energetic minimax of the paths of the Kawasaki dynamics from 
□ to ffl, which all have energy T* because H{D) = 0. 

In [3] we showed that the first entrance distribution on the set of critical droplets is uniform, 
computed the expected transition time up to and including a multiplicative factor of order one, and 
proved that the nucleation time divided by its expectation is exponentially distributed, all in the 
limit as /3 — >■ cx3. These results, which are typical for metastable behavior, were proved under three 
hypotheses: 

(HI) A'stab = ffl- 

(H2) There exists aV* kT" such that < V* for all 77 e X\{a, ffl}. 

(H3) A hypothesis about the shape of the configurations in the essential gate for the transition from 
□ to ffl (for details see [3]). 

Hypotheses (H1-H3) are the geometric input that is needed to derive the main theorems in pj] with 
the help of the potential-theoretic approach to metastability as outlined in Bovier [2J. In the present 
paper we prove (II1-II2) and identify the energy of critical droplets. In [4J we settle the rest. 
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Lemma 1.2 (H1-H2) imply that V\j = T* , and hence that A'^ota = D ind T = T* . 

Proof. By Definition [TTIe-h) and (HI), ffl G which imphes that Vq < T*. We show that (H2) 
implies Vu = F*. The proof is by contradiction. Suppose that Vb < F*. Then, by Definition ll.lf h). 
there exists a rjQ £ In\ffl such that ^(□,7yo) — H{n) < T*. But (H2), together with the finiteness of 
X , imphes that there exist an m G N and a sequence 771, ... , rjm € X with rjm = ffl such that ryj+i G I,,^ 
and <?(77i, ryi+i) < H{rii) + V* for i = 0, . . . , m — 1. Therefore 

<?(?7o, ffl) < max ^>(?7„ r/,+i) < max [H{f^,) + V*] = H{r]o) + ¥" < H{n) + T\ (1.9) 

t— 0,...,m— 1 i— 0,...,m— 1 

where in the first inequahty we use that <I>{ri,a) < max{<?(77, ^), <?(^, ct)} for aU ri,a,^ G X, and in the 
last inequality that rjQ E and V* < F*. It follows that 

<?(□, ffl) - H{D) < max{<Z'(n, 770) - H{n), <P{r]a, ffl) - H{D)} < F*, (1.10) 

which contradicts Definition 1 1 . 1 f h) . Observe that the proof uses that ^'mota consists of a single config- 
uration. ■ 

Hypotheses (H1-H2) imply that (A'meta, '%stab) = (n,ffl)) and that the highest energy barrier be- 
tween any two configurations in X is the one separating □ and ffl, i.e., (□, ffl) is the unique metastahle 
pair. Hypothesis (H3) is needed only to find the asymptotics of the prefactor of the expected transition 
time in the limit as A — > Z^. The main theorems in [3] involve three model- dependent quantities: the 
energy, the shape and the number of critical droplets. 



1.3 Main theorems 

In [3] it was shown that Ai + A2 < 4[/ is the metastahle region, i.e., the region of parameters for 
which □ is a local minimum but not a global minimum of H. Moreover, it was argued that within this 
region the subregion where Ai, A2 < U is of no interest because the critical droplet consists of two 
free particles, one of type 1 and one of type 2. Therefore the proper metastahle region is 



< Ai < A2, Ai + A2 < 4J7, A2 > [/, 



fl.ll) 



as indicated in Fig. [T] 




Figure 1: Proper metastahle region. 



In this present paper, the analysis will be carried out for the subregion where 

0<Ai<C/, A2-Ai>2[/, Ai+A2<4C/, (1.12) 

as indicated in Fig. [5J Note: The second and third restriction imply the first restriction. Nevertheless, 
we write all three because each plays an important role in the sequel. 

The following three theorems are the main result of the present paper and are valid subject to 
p.l2p . We write [•] to denote the upper integer part. 
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Az-Ai =2U 



A, 

Figure 2: Subregion of the proper metastable region given by (jl.l2p . 
Theorem 1.3 Xsti,h = H- 

Theorem 1.4 There exists aV* < lOU — Ai such that Vjj < V* for all rj e A'\{n, ffl}. Consequently, 
ifV* > IOC/ - Ai, then Af^ota = □ andV ^T*. 



Theorem 1.5 T* = -[£*{£* - 1) + 1](4C/ - Aj - Aa) + (2£* + l)Ai + Aa wit/i 

Ai 



4C/- Ai - As 



e N. 



(1.13) 



Theorem 11.31 settles hypothesis (HI) in i3j, Theorem II .41 settles hypothesis (H2) in ^3^ when T* > 
lOU - Ai, while Theorem O identifies T*. 

As soon as V* < T*, the energy landscape does not contain wells deeper than those surrounding □ 
and ffl. Theorems 11.31 and 11.41 imply that this occurs at least when F* > lOU ~ Ai, while Theorem II .51 
identifies F* and allows us to exhibit a further subregion of (jl.l2|) where the latter inequality is satisfied. 
This further subregion contains the shaded region in Fig. |31 




Ai 

Figure 3: The parameter region where F* > lOU — Ai contains the shaded region. 



1.4 Discussion 

1. In Section 21 we will see that the critical droplets for the crossover from □ to ffl consist of a rhombus- 
shaped checkerboard with a protuberance plus a free particle, as indicated in Fig. SI A more detailed 
description will be given in [4j. 

2. Abbreviate 

e = 4t/-Ai-A2 (1.14) 
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Figure 4: A critical droplet. Light-shaded squares are particles of type 1, dark-shaded squares are 
particles of type 2. The particles of type 2 form an 1* x {£* — 1) quasi-square with a protuberance 
attached to one of its longest sides, and are all surrounded by particles of type 1. In addition, there is 
a free particle of type 2. As soon as this free particle attaches itself "properly" to a particle of type 1 
the dynamics is "over the hill" (see [3], Section 2.3, item 3). 

and write = (Ai/e) -|- l with t g [0, 1). Then an easy computation shows that T* — (Ai)^/e + Ai + 
411 + £i(l — i). From this we see that 

r~Ai/e, r*-(Ai)Ve, £ 10- (1-15) 

The limit e J, corresponds to the weakly supersaturated regime, where the lattice gas wants to 
condensate but the energetic threshold to do so is high (because the critical droplet is large) . From the 
viewpoint of metastability this regime is the most interesting. The shaded region in Fig. [3] captures 
this regime for all < Ai < [/. This region contains the set of parameters where (Ai)^/£ + Ai + AU > 
lOU - Ai, i.e., e/U < (Ai/C/)V[6 - 2(Ai/C/)]. 

3. The simplifying features of (|1.12p over (|l.lip are the following: Ai < J7 implies that each time a 
particle of type 1 enters A and attaches itself to a particle of type 2 in a droplet the energy goes down, 
while A2 — Ai > 2U implies that no particle of type 2 sits on the boundary of a droplet that has 
minimal energy given the number of particles of type 2 in the droplet. In [3] we conjectured that the 
metastability results presented there actually hold throughout the region given by (jl.lip , even though 
the critical droplets will be different when Ai > U. 

As will become clear in Section [31 the constraint Ai < U has the effect that in all configurations 
that are local minima of H all particles on the boundary of a droplet are of type 1. It will turn out 
that such configurations consist of a single rhombus-shaped checkerboard droplet. We expect that as 
Ai increases from U to 2U there is a gradual transition from a rhombus-shaped checkerboard critical 
droplet to a square-shaped checkerboard critical droplet. This is one of the reasons why it is difficult 
to go beyond (ITTT^ . 

4. What makes Theorem II .41 hard to prove is that the estimate on has to be uniform in rj ^ {□, ffl}. 
In configurations containing several droplets and/or droplets close to d~A there may be a lack of free 
space making the motion of particles inside A difficult. The mechanisms developed in Section [5] allow us 
to realize an energy reduction to a configuration that lies on a suitable reference path for the nucleation 
within an energy barrier 10?7 — Ai also in the absence of free space around each droplet. 

We will see in Section [5] that for droplets sufficiently far away from other droplets and from 9^ A 
a reduction within an energy barrier < AU + Ai is possible. Thus, if we would be able to control the 
configurations that fail to have this property, then we would have V* < AU + Ai and, consequently, 
would have A'mcta — D and F = F* throughout the subregion given by (jl.l2p because F"^ > AU -I- Ai. 

Another way of phrasing the last observation is the following. We view the "liquid phase" as the 
configuration filling the entire box A. If, instead, we would let the liquid phase correspond to the set 
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of configurations filling most of A but staying away from d A, then the metastability results derived 
in [3] would apply throughout the subregion given by (jl.l2p . 

5. Theorems 11.31 and 11.51 can actually be proved without the restriction A2 — Ai > 2U . However, 
removal of this restriction makes the task of showing that in droplets with minimal energy all particles 
of type 2 are surrounded by particles of type 1 more involved than what is done in Section [31 We omit 
this extension, since the restriction A2 — Ai > 2U is needed for Theorem II .41 anvwav. 

Outline. Section [2] contains preparations. Theorems 1 1 . BHTTSl are proved in Sections [3HS1 respectively. 
The proofs are purely combinatorial^ and are rather involved due to the presence of two types of 
particles rather than one. Sections [SHU deal with statics and Section [S] with dynamics. Section [5] is 
technically the hardest and takes up about half of the paper. More detailed outlines are given at the 
beginning of each section. 



2 Coordinates, definitions and polyominoes 

Section 12.11 introduces two coordinate systems that are used to describe the particle configurations: 
standard and dual. Section 12.21 lists the main geometric definitions that are needed in the rest of the 
paper. Section [2.31 proves a lemma about polyominoes (finite unions of unit squares) and Section [2.4l a 
lemma about 2-tiled clusters (checkerboard configurations where all particles of type 2 are surrounded 
by particles of type 1). These lemmas are needed in Section [3] to identify the droplets of minimal 
energy given the number of particles of type 2 in A. 



2.1 Coordinates 

1. A site j e A is identified by its standard coordinates [xi{i), X2{i)), and is called odd when xi(i)+X2(i) 
is odd and even when xi{i) + x^ii) is even. The standard coordinates of a particle p in A are denoted 
by x(p) = {xi{p), X2{p))- The parity of a particle p is defined as xi{p) + X2{p) + r]{x{p)) modulo 2, and 
p is said to be odd when the parity is 1 and even when the parity is 0. 

2. A site z € A is also identified by its dual coordinates 

xi{i)-X2ii) ... xi{i)+X2{i) , . 

uiW = 2 ' ^ ^ 2 ' ^ ' 

Two sites i and j are said to be adjacent, written i ^ j, when |a::i(«) — xi{j)\ + |a;2(«) — X2{j)\ = 1 or, 
equivalently, \ui{i) - ui{j)\ = \u2{i) - U2{j)\ = | (see Fig. [5]). 

3. For convenience, we take A to be the (L + |) x (L + |) dual square centered at the origin for some 
L e N with L > 2£* (to allow for H(S) < H{n); see Section [3T|) . Particles interact only inside A", 
which is the {L + x {L + ^) dual square centered at the origin. This dual square, a rhombus in 
standard coordinates, is convenient because the local minima of H are rhombus-shaped as well (see 
Section [3]). 



2.2 Definitions 

1. A site « e A is said to be lattice- connecting in the configuration 77 if there exists a lattice path A 
from i to d~A such that 77 (j) = for all j £ X with j i. We say that a particle p is lattice-connecting 
if x{p) is a lattice-connecting site. 

2. Two particles in 77 at sites i and j are called connected if i ^ j and ri(i)rj{j) = 2. If two particles pi 
and p2 are connected, then we say that there is an active bond b between them. The bond b is said to 
be incident to pi and p2. A particle p is said to be saturated if it is connected to four other particles, 
i.e., there are four active bonds incident to p. The support of the configuration 77, i.e., the union of the 
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(a) (b) 



Figure 5: A configuration represented in: (a) standard coordinates; (b) dual coordinates. Light-shaded 
squares are particles of type 1, dark-shaded squares are particles of type 2. In dual coordinates, particles 
of type 2 are represented by larger squares than particles of type 1 to exhibit the "tiled structure" of 
the configuration. 



unit squares centered at the occupied sites of 77, is denoted by supp(77). For a configuration 77, ^.1(77) 
and 712(77) denote the number of particles of type 1 and 2 in 77, and i?(ry) denotes the number of active 
bonds. The energy of 77 equals H{ri) = Aini{ri) + A2'n2(r]) — UB(ri). 

3. Let G{ri) be the graph associated with j], i.e., G{ri) = {V{r]), E{r])), where V{ri) is the set of sites 
i G A such that 77(2) 7^ 0, and E{7]) is the set of the pairs i,j G V{ri), such that the particles 
at sites i and j are connected. A configuration ?/ is called a subconfiguration of rj, written rj' -< rj, if 
r]'{i) = ri{i) for all i G A such that r]'{i) > 0. A subconfiguration c ^ 77 is a cluster if the graph G(c) is 
a maximal connected component of G{ri). The set of non-saturated particles in c is called the boundary 
of c, and is denoted by dc. Clearly, all particles in the same cluster have the same parity. Therefore 
the concept of parity extends from particles to clusters. 

4. For a site i e A, the tile centered at i, denoted by t{i), is the set of five sites consisting of i and the 
four sites adjacent to i. If i is an even site, then the tile is said to be even, otherwise the tile is said 
to be odd. The five sites of a tile are labeled a, b, c, d, e as in Fig. [HI The sites labeled a, b, c, d are 
called junction sites. If a particle p sits at site i, then t{i) is also denoted by t(p) and is called the tile 
associated with p. In standard coordinates, a tile is a square of size v^- In dual coordinates, it is a 
unit square. 

5. A tile whose central site is occupied by a particle of type 2 and whose junction sites are occupied by 
particles of type 1 is called a 2-tile (see Fig. [5]). Two 2-tiles are said to be adjacent if their particles of 
type 2 have dual distance 1. A horizontal (vertical) 12-bar is a maximal sequence of adjacent 2-tiles 
all having the same horizontal (vertical) coordinate. If the sequence has length 1, then the 12-bar is 
called a 2-tiled protuberance. A cluster containing at least one particle of type 2 such that all particles 
of type 2 are saturated is said to be 2-tiled. A 2-tiled configuration is a configuration consisting of 
2-tiled clusters only. 






■ 





-0 





■ ■ ■ 










■ 





-0 ■- 


(a) 


(b) 


(c) 





(d) 



Figure 6: Tiles: (a) standard representation of the labels of a tile; (b) standard representation of a 
2-tile; (c) dual representation of the labels of a tile; (d) dual representation of a 2-tile. 



8 



6. The tile support of a configuration 77 is defined as 

W= U (2.2) 

pero2(»;) 

where 1372 is the set of particles of type 2 in ?/. Obviously, [77] is the union of the tile supports of 
the clusters making up 77. For a standard cluster c the dual perimeter, denoted by P(c), is the length 
of the Euclidean boundary of its tile support [c] (which includes an inner boundary when c contains 
holes). The dual perimeter P{ri) of a 2-tiled configuration 77 is the sum of the dual perimeters of the 
clusters making up 77. 

7. V*^„2 is the set of configurations such that in A the number of particles of type 2 is 712. V^"^^ 
the set of configurations such that in A the number of particles of type 2 is 712, the number of active 
bonds is 4^2 , and there is no isolated particle of type 1. In other words, V^^^ is the set of 2-tiled 
configurations with 77,2 particles of type 2. The lower index -k is used to indicate that configurations in 
these sets can have an arbitrary number of particles of type f . A configuration 77 is called standard if 
77 G V^"^^, and its tile support is a standard polyomino in dual coordinates (see Definition 12. f I below 
for the definition of a standard polyomino). 

8. A unit hole is an empty site such that all four of its neighbors are occupied by particles of the same 
type (either all of type f or all of type 2). An empty site with three neighboring sites occupied by 
a particle of type f is called a good dual corner. In the dual representation a good dual corner is a 
concave corner (see Fig. [7]). 

2.3 A lemma on polyominoes 

The tile support of a cluster c can be represented by a polyomino, i.e., a finite union of unit squares. 
The following notation is used: 

£i{c) = width of c (= number of columns). 
£2(0) = height of c (= number of rows). 

Vi (c) = number of vertical edges in the z-th non-empty row of c. 

hj (c) = number of horizontal edges in the j-th non-empty column of c. 

P(c) = length of the perimeter of c. 

(5(c) — number of holes in c. 

ipic) — number of convex corners of c. 

0(c) = number of concave corners of c. 

Note that ip^c) — '^i'Ji^ '/'(^) ^^nd 0(c) — 'J2i!^i^ 0(*)j where A^(c) is the number of vertices in the 
polyomino representing c. If two edges ci and 62 are incident to vertex i at a right angle with a unit 
square inside and no unit squares outside, then •ip{i) = 1 and = (Fig. [71^ a)). On the other hand, 
if there is no unit square inside and three unit squares outside, then 7/1(1) = and ^(i) = I (Fig.jTJb)). 
If four edges ci, e2, 63, 64 are incident to vertex i, with two unit squares in opposite angles, then 
Vj(i) = and 0(i) = 2 (Fig. [Tjc)). 

Definition 2.1 [Alonso and Cerf pj.] A polyomino is called monotone if its perimeter is equal to 
the perimeter of its circumscribing rectangle. A polyomino whose support is a quasi-square (i.e., a 
rectangle whose side lengths differ by at most one), with possibly a bar attached to one of its longest 
sides, is called a standard polyomino. 
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(a) (b) (c) 

Figure 7: Corners of polyominoes: (a) one convex corner; (b) one concave corner; (c) two concave 
corners. Shaded mean occupied by a unit square. 

In the sequel, a key role will be played by the quantity 

r(c) = 2P(c) + [V'(c) - (f>{c)] = 2P(c) + 4 - 4Q(c). (2.3) 

Lemma 2.2 (i) All polyominoes c with a fixed number of monominoes minimizing T{c) are single- 
component monotone polyominoes of minimal perimeter, which include the standard polyominoes. 
(a) If the number of monominoes is P' , — \, £{£ — 1) or £{i — 1) — 1 for some £ G N\{1}, then the 
standard polyominoes are the only minimizers ofT{c). 

Proof. In the proof we assume w.l.o.g. that the polyomino consists of a single cluster c. 

(i) The proof uses projection. Pick any non-monotone cluster c. Let 

£= (^2 0^i)(c), (2.4) 

where TT2 and tti denote the vertical, respectively, the horizontal projection of c. The effect of vertical 
and horizontal projection is illustrated in Fig. [51 By construction, c is a monotone polyomino (see e.g. 
the statement on Ferrers diagrams in the proof of Alonso and Cerf [T], Theorem 2.2). 
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Figure 8: Effect of vertical and horizontal projection. 

Suppose first that Q{c) = 0. Then T(c) = 2P(c)-t-4. Since c is not monotone, we have P{c) < P{c), 
and so c is not a minimizer of T(c) . 

Suppose next that Q{c) > 1. Since 

l2{c) hie) 

P(c)= ^x;,(c)+^/i,(c) (2.5) 

-i=i j=i 

and every hole belongs to at least one row and one column, we have 

P(c) > 2[4(c) + £2{c)] + 4Q(c). (2.6) 
On the other hand, since c is a monotone polyomino, we have Vi{c) = hj{c) = 2 for all i and j, and so 

P(5) = 2[£i(2)+^2(5)]. (2.7) 
Moreover, since £i{c) < £i{c) and £2(2) < £2(0), we can combine (j2.6H2.7l) to get 

P{£)-P{c)<-AQ{c), (2.8) 
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Using (j2.8p . we obtain 



r{£) - r{c) = [2P{c) + 4] - [2P(c) + 4 - 4g(c)] = 2[P(5) - P(c)] + 4Q(c) < -4Q(c) < -4 < 0, (2.9) 
and so c is not a minimizcr of T(c). 

(ii) We saw in the proof of (i) that if c is a minimizer of T(c), then c is monotone, and hence does not 
contain holes and minimizes P{c). The claim therefore follows from Alonso and Cerf [Ij, Corollary 3.7, 
which states that if the number of monominoes is — 1, i{£ — 1) or £{£— 1) — 1 for some £ E N\{1}, 
then the standard polyominoes are the only minimizers of P{c). I 

2.4 Relation between T and the number of missing bonds in 2— tiled clusters 

In this section we consider 2-tiled clusters and link the number of particles of type 1 and type 2 to the 
number of active bonds and the geometric quantity T considered in Section 12.31 

Lemma 2.3 For any 2-tiled cluster c (i.e., c G V^"^ for some n2), 4ni(c) — B{c)+T{c) and 4n2{c) — 
B{c). 

Proof. The claim of the lemma is equivalent to the affirmation that T(c) = M{c) with M{c) the 
number of missing bonds in c. Indeed, informally, for every unit perimeter two bonds are lost with 
respect to the four bonds that would be incident to each particle of type 1 if it were saturated, while 
one bond is lost at each convex corner and one bond is gained at each concave corner. 

Formally, let p be a particle of type 1, B{p) the number of bonds incident to p, and M [p) = A— B{p) 
the number of missing bonds of p. Consider the set of particles of type 1 at the boundary of a 2-tiled 
cluster, i.e., the set of non-saturated particles of type 1. Each of these particles belongs to one of four 
classes (see Fig. IH]): 

class 1: p has two neighboring particles of type 2 belonging to the same 12-bar. 
class 2: p has two neighboring particles of type 2 belonging to different 12-bars. 
class 3: p has three neighboring particles of type 2. 
class 4: p has one neighboring particle of type 2. 




(a) (b) (c) (d) 



Figure 9: The circled boundary particle of type 1 belongs to: (a) class 1; (b) class 2; (c) class 3; (d) 
class 4. 

Let Mk (c) be the number of missing bonds of particles of class k in cluster c, and Ak (c) the number 
of edges incident to particles of class k in cluster c. Then 

Mi(c) = 2, Ai(c) =2; Af2(c) = 2, A2(c)=4; A-h{c) ^ I , A^{c) = 2- M4(c) = 3, ^4(0) = 2. 

(2.10) 

Let Nk{c) be the number of particles of class k of type 1 in cluster c. Observing that a cluster has 
two concave corners per particle of class 2. one concave corner per particle of class 3 and one convex 
corner per particle of class 4, we can write 

r(c) = 2P{c) - 2N2{c) - Nsic) + mic). (2.11) 
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Since the dual perimeter of a cluster is equal to its total number of dual edges, we have 

4 

2P{c)^Y.Ak{c)Nk{c) = 2Ni{c)+m2{c)+2N3{c) + 2Ni{c) (2.12) 

k=l 

(the sum counts each edge of the 2-tile twice). The total number of missing bonds, on the other hand, 
is 

4 

M(c) = Mk{c)Nk{c) = 2iVi(c) + 2iV2(c) + N^ic) + SN^ic). (2.13) 

k=l 

Combining (|2JTH2?T3)l . we arrive at T{c) M{c). ■ 

3 Proof of Theorem II. 3t identification of A^tab 

Recall that (the part of A where particles interact) is an (i + |) x (L + i) dual square with 
L > 2£*. Let 77stab, f^stab Configurations consisting of a 2-tiled dual square of size L with even 

parity, respectively, odd parity. These two configurations have the same energy. Theorem II .31 says that 
■^stab = {^stab, ?7stab} ~ ^- Scctiou [3. II Contains two lemmas about 2-tiled configurations with minimal 
energy. Section uses these two lemmas to prove Theorem 11.31 

3.1 Standard configurations are minimizers among 2— tiled configurations 

Lemma 3.1 Within Vf"^^^; standard configurations achieve the minimal energy. 

Proof. Recall from item [5] in Section that 

Hi?^) = Aini(7?) + A2n2{v) - UB{r,). (3.1) 

In V^'Jj^2 both and B — 4?i2 are fixed, and hence min^gy4„2 H{rj) is attained at a configuration 
mmimizmg ni. By Lemma [131 if 77 e V^%, then 

ni{ri) = \[B{r,) + T{v)], min) - \B{fj). (3.2) 

Hence, to minimize ni{r]) we must minimize T{vi). The claim therefore follows from Lemma l2.2f i). ■ 

For a standard configuration the computation of the energy is straightforward. For ^ G N, C G {0, 1} 
and fc G No with fc < f + C, let rj^^^^^ denote the standard configuration consisting of an x + C,) 
(quasi-)square with a bar of length k attached to one of its longest sides (see Fig. [TUl) . 



Figure 10: A standard configuration with £ = 7, ^ = 1 and k = 5. 

Lemma 3.2 The energy of rj^'^'^ is (recall p.l4p ) 

Hij^''^^'') = -eW + C) + fc] + Ai + (^ + C) + 1 + l{fc>o}]. (3.3) 
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Proof. Note that P{f^^-^^^) = 2\l + + C) + l{fe>o}] and Qi.jf-'^-^) = 0, so that 

rCry^^'^'^-) = 4[£ + + C) + 1 + l{fc>o}]- (3.4) 

Also note that 

B(r/^'^''=) = 4[£+(^ + C) + A:], (3.5) 
because ah particles of type 2 are saturated. However, by (l3.1H3.2p . we have 

7j(^^,C,fe) ^ -\eB{^'^^■^) + ir(r7^-^''=)Ai, (3.6) 

and so the claim follows by combining p.4H3.6p . ■ 

Note that the energy increases by Ai — e (which is > if and only if > 2 by (|1.13p l when a bar 
of length A: = 1 is added, and decreases by e each time the bar is extended. Note further that 

7?(r;^^i'0) - ^ _ i7(r;^+i^0'0) - i7(r;^'i'0) ^ Ai - (^ + l)e, (3.7) 

which show that the energy of a growing sequence of standard configurations goes up when i < 1* and 
goes down when i > . The highest energy is attained at if "^'^^i, which is the critical droplet in 

Fig.m 

It is worth noting that H{ri^^ '''''^) < 0, i.e., the energy of a dual square of side length 2£* is lower 
than the energy of □. This is why we assumed L > 2£*, to allow for -ff(ffl) < H(D). 



3.2 Stable configurations 

In this section we use Lemmas I3.1H3.2I to prove Theorem 11.31 

Proof. Let rj denote any configuration in Xstah- Below we will show that: 

(A) rj does not contain any particle in A. 

(B) 77 is a 2-tiled configuration, i.e., 77 G V^'^^^ some 712 (= n2{rj)). 

Once we have (A) and (B), we observe that 77 cannot contain a number of 2-tiles larger than L^. 
Indeed, consider the tile support of 77. Since A~ is an {L + ^) x (L + i) dual square, if the tile support 
of 77 fits inside A~, then so does the dual circumscribing rectangle of rj. But any rectangle of area > 
has at least one side of length L + 1. Hence 71.2(7;) < L^, and therefore the number of 2-tiles in 77 is 
at most L^. By Lemmas I3.1H3.21 the global minimum of the energy is attained at the largest dual 
quasi-square that fits inside A~, since L > 2i*. We therefore conclude that 77 G {??stab, ??stab}' which 
proves the claim. 

Proof of (A). Since in A particles do not feel any interaction but have a positive energy cost, removal 
of a particle from A always lowers the energy. 

Proof of (B). We note the following three facts: 

(1) rj does not contain isolated particles of type 1. 

(2) d^A~ does not contain any particle of type 2. 

(3) All particles of type 2 in 7; have all their neighboring sites occupied by a particle. 

For (1), simply note that the configuration obtained from 77 by removing isolated particles has lower 
energy. For (2), note that particles in d~A~ have at most two active bonds. Therefore, if 77 would 
have a particle of type 2 in d~A~, then the removal of that particle would lower the energy, because 
A2 - Ai > 2U and Ai > (recall (fTTT^ 'l imply A2 > 2U. For (3), note that if a particle of type 2 
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has an empty neighboring site, then the addition of a particle of type 1 at this site lowers the energy, 
because Ai < U (recall (|1.12p l. 

We can now complete the proof of (B) as follows. The constraint A2 — Ai > 2U implies that 
any particle of type 2 in 77 must have at least three neighboring sites occupied by a particle of type 
1. Indeed, the removal of a particle of type 2 with at most two active bonds lowers the energy. 
But the fourth neighboring site must also be occupied by a particle of type 1. Indeed, suppose that 
this site would be occupied by a particle of type 2. Then this particle would have at most three 
active bonds. Consider the configuration 77 obtained from rj after replacing this particle by a particle 
of type 1. Then B{fj) — B{rj) > —2, 71,1(77) — 7^1(7;) — 1 and 712(77) — 77,2(77) = —1. Consequently, 
H{rj) — H{ri) < Ai — A2 + 2U < 0. Hence, any particle of type 2 in 77 must be saturated. ■ 

4 Proof of Theorem 11.51 : identification of = S) 

In Section HTTI we prove Theorem 11.51 sub iect to the following lemma. 

Lemma 4.1 For any 77,2 < , the configurations of minimal energy with 772 particles of type 2 belong 
to V^"j^^, i.e., are 2-tilcd configurations. 

The proof of this lemma is given in Section [4.21 

4.1 Proof of Theorem 11.51 subject to Lemma 14.11 

Proof. For y C X, define the external boundary of y by dy = {77 G X\y: 3ri' E y, rj ^ rj'} and 
the bottom of y by J-{y) — argmin^gj; H{rj). According to Manzo, Nardi, Olivieri and Scoppola [5J, 
Section 4.2, ^(D, ffl) = min^ggg H{rf) ioi B C X any (!) set with the following properties: 

(I) B is connected via allowed moves, D E B and Si ^ B. 
(II) There is a path w* : □ ^- ffl such that {argmax,,g^* H{r])} n TidB) ^ 0. 

Thus, our task is to find such a B and compute the lowest energy of dB. 

For (I), choose B to be the set of all configurations 77 such that n2{ri) < ^*(^* — 1) + 1. Clearly this 
set is connected, contains □ and does not contain ffl. 

For (II), choose lu* as follows. A particle of type 2 is brought inside A {AH = A2), moved to the 
origin and is saturated by four times bringing a particle of type 1 (AH = Ai) and attaching it to the 
particle of type 2 (AH = —U). After this first 2-tile has been completed, follows a sequence of 
increasing 2-tiled dual quasi-squares. The passage from one quasi-square to the next is obtained by 
adding a 12-bar to one of the longest sides, as follows. First a particle of type 2 is brought inside A 
(AH — A2) and is attached to one of the longest sides of the quasi-square (A_ff = —2U). Next, twice 
a particle of type 1 is brought inside the box {AH = Ai) and is attached to the (not yet saturated) 
particle of type 2 {AH — —U) in order to complete a 2-tiled protuberance. Finally, the 12-bar 
is completed by bringing a particle of type 2 inside A {AH = A2), moving it to a concave corner 
{AH = — 3[/), and saturating it with a particle of type 1 {AH = Ai, respectively, AH = —U). It is 
obvious that uj* eventually hits ffl. The path uj* is referred to as the reference path for the nucleation. 

Call 7/* the configuration in uj* consisting of an (* y. (^* — 1) quasi-square, a 2-tiled protuberance 
attached to one of its longest sides, and a free particle of type 2 (see Fig. [TTJ there are many choices 
for UJ* depending on where the 2-tiled protuberances are added; all these choices are equivalent. Note 
that, in the notation of Lemma [3.21 77* — rj^ -i>i'i -|- fp[2], where -Hfp[2] denotes the addition of a free 
particle of type 2. Observe that: 

(a) UJ* exits B via the configuration rj*; 

(b) 77* e F{dB); 
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(c) 77* e {argmax^g,^* H{i])}. 



Observation (a) is obvious, while (b) follows from Lemmas 13.11 and 14.11 To see (c), note the following: 
(1) The total energy difference obtained by adding a 12-bar of length £ on the side of a 2-tiled cluster 
is Ai7(adding a 12-bar) = Ai — e£, which changes sign at ^ = ^* (recall (|3.7p ); (2) The configurations 
of maximal energy in a sequence of growing quasi-squares are those where a free particle of type 
2 enters the box after the 2-tiled protuberance has been completed. Thus, within energy barrier 
2Ai + 2A2 — 4C/ = 4[/ — e the 12-bar is completed downwards in energy. This means that, after 
configuration 77* is hit, the dynamics can reach the 2-tiled dual square of £* x £* while staying below 
the energy level H{r]*). Since all 2-tiled dual quasi-squares larger than £* y. {£* — 1) have an energy 
smaller than that of the 2-tiled dual quasi-square £* y. {£* — \) itself, the path does not again reach 
the energy level H(rj*). 

Because of (a-c), we have ffl) = H{'ri*). To complete the proof, use Lemma to compute 



-14,1 



+ fp[2]) = -e[r(r - 1) + 1] + Ai(2r + 1) + A2 



(4.1) 



Figure 11: A critical configuration rj* . This is the dual version of the critical droplet in Fig. 



4.2 Proof of Lemma 14.11 

The proof of Lemma HTTI is carried out in two steps. In Section 14.2.1 1 we show that the claim holds for 
single-cluster configurations with a fixed number of particles of type 2. In Section 14.2.21 we extend the 
claim to general configurations with a fixed number of particles of type 2. 

4.2.1 Single clusters of minimal energy are 2 tiled clusters 

Lemma 4.2 For any single-cluster configuration rj £ '^*,n2\'^*"n2 there exists a configuration fj £ V^"^^ 
such that H{fj) < H^rj). 

Proof. Pick any r/ e V*_„2\Vf'"^2. Every neighboring site of a particle of type 2 in the cluster is either 
empty or occupied by a particle of type 1, and there is at least one non-saturated particle of type 2. 
Since 77 consists of a single cluster, fj can be constructed in the following way: 

• fjii) — rj(i) for all i £ supp(77). 

• vii) — 1 for all j (f. supp(77) such that there exists an i ^ j with 77(7) — 2. 
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Since 

H{ri) = Aini(r/) + AanaCr,) - UB{ri), 
H{f,) = Aini(r)) + ^^n^iv) - UB{f,), 

and n2(ji) = 71.2(77), we have 

H(f,) - i/(77) = Ai[77i(7y) - 771(77)] - U[B(f,) - ^(77)]. (4.3) 

By construction, ^(77) - ^(77) > 77.1(77) - rti(7;) > 0. Since < Ai < U (recall (|1.12p '). it follows from 
(itill) that H{fi) < H{ri). ■ 



4.2.2 Configurations of minimal energy with fixed number of particles of type 2 

Lemma 4.3 For any 712 and any configuration rj G V*_„2 consisting 0/ at least two clusters, any 
configuration rf such that rj* is a single cluster, rf £ Vfn^ '^^'^ V* ^ standard configuration satisfies 
H{rf)<H{r]). 

Proof. Let 77 e V*,n2 be a configuration consisting of /c > 1 clusters, labeled ci, . . . ,Ck- Let rj"'^'^'^^^ 
denote any standard configuration with 712(0^) particles of type 2. By Lemmas 13.11 and 14.21 we have 

H(77) = ^i/(c,)>^W^(^-)). (4.4) 

i=l 1=1 

By Lemma [2.31 we have (recall (jl.l4p ) 

k k 

i/(77"=(=-)) = [Ami (77"^ ("•^) + A2772(77"^(^')) - ?7B(7;"^(=-))] 

i=l i=l 
fe 

= ^ [Al{7l2(77"^(=')) + ir(77"^(^-))} + A2n2(77"-(=')) - f/47l2(77"^(='^)] (4.5) 
i=l 

fc 

= ^ [ - £712(77"^("'^) + iAir(77"^("-))] . 

1=1 

But from Lemma [52] it follows that 

fc 

i=l 

where "2('=i) denotes any standard configuration with X]i=i"-2(ci) = 712(77) particles of type 2. 

Combining (|4.4H4.6p . we arrive at 

H{v) > -£712(77) + iAir(77"=('')) = ff(77"^('')). (4.7) 



5 Proof of Theorem II. 4t upper bound on for 77 ^ {□, ffl} 

In this section we show that for any configuration rj ^ {□, ffl} it is possible to find a path ui: rj ^ rj' 
with 7y' e {□,ffl} such that max^g^ ff(^) < H{ri) + V* with V* < lOU - Ai and -q' e By 
Definition 11.1 f c-e) . this implies that < V* for all 77 ^ {□, ffl} and therefore settles Theorem II .41 

Section TS .31 describes an energy reduction algorithm to find uj. Roughly, the idea is that if rj contains 
only "subcritical clusters", then these clusters can be removed one by one to reach □, while if 7/ contains 
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some "supercritical cluster", then this cluster can be taken as a stepping stone to construct a path to 
ffl that goes via a sequence of increasing rectangles. In particular, the supercritical cluster is first 
extended to a 2-tiled rectangle touching the north-boundary of A, after that it is extended to a 2-tiled 
rectangle touching the west-boundary and the east-boundary of A, and finally it is extended to ffl. 

To carry out this task, six energy reduction mechanisms are needed, which are introduced and 
explained in Section 15.21 

• Moving unit holes inside 2-tiled clusters (Section 15.2. ip . 

• Adding and removing 12-bars from lattice-connecting rectangles (Section 15. 2. 2p . 

• Changing bridges into 12-bars (Section 15. 2. 3p . 

• Maximally expanding 2-tiled rectangles (Section 15. 2. 41) . 

• Merging adjacent 2-tiled rectangles (Section 15. 2. 5[) . 

• Removing subcritical clusters (Section 15. 2. 6p . 

Each of Sections I5.2.1H5^2.6I states a definition and a lemma, and uses these to prove a proposition 
about the relevant energy reduction mechanism. The six propositions thus obtained will be crucial for 
the energy reduction algorithm in Section [5.31 

In Section [5. II we begin by defining beams and pillars, which are needed throughout Section [5.21 
5.1 Beams and pillars 

Lemma 5.1 Let r/ be a configuration containing a tile t that has at least three junction sites occupied 
by a particle of type 1. Then the configuration rj' obtained from by turning t into a 2-tile satisfies 
H{rj')<H{r^). 




(i) (ii) (iii) (iv) (v) (vi) (vii) (viii) 



Figure 12: Possible tiles with at least three junction sites occupied by a particle of type 1. 

Proof. W.l.o.g. we may assume that ri{ta) = rj{ti,) = rj{td) = 1, and that r]' is the configuration in 
Fig. jSJd), i.e., rj'{ta) = 'n'iib) = v'i^c) = v'{td) = 1, 77'(^e) = 2. The following eight cases are possible 
(see Fig. inland recaU pTT^ I: 

(i) {ri{tc),ri{te)) = (0,0). One particle of type 1 and one particle of type 2 are added, and at least 
four new bonds are activated: AH < Ai -I- A2 — 4C7 < 0. 

(ii) {ri{tc),ri{te)) — (0,2). One particle of type 1 is added, and one new bond is activated: AH = 
Ai-U <0. 

(iii) {ri{tc) , ri{t(,)) — (2,0). One particle of type 2 is moved to another site without deactivating any 
bonds, after which case ^ applies. 

(iv) {ri{tc),ri{te)) — (2,2). One particle of type 2 with at most three active bonds is replaced by one 
particle of type 1 with at least one active bond: AH < Ai — A2 + 2U < 0. 

(v) {ri{tc),ri{te)) — (1,0). One particle of type 2 is added, and four new bonds are activated: AH = 
A2 - 4J7 < 0. 
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(vi) {ii{tc),ri{te)) — (0,1). One particle of type 1 is moved to another site without deactivating 
any active bond, one particle of type 2 is added, and at least four new bonds are activated: 
AiJ < A2 - 4[/ < 0. 

(vii) {ri{tc),ri{te)) — (2, 1). Two particles are exchanged without deactivating any bonds: AH < 0. 

(viii) {rj{tc),rj{te)) = (1,1)- One particle of type 1 is replaced by a particle of type 2, and four new 
bonds are activated: /S.H = A2 — Ai — AU < 0. 



Definition 5.2 A beam of length I is a row (or column) of £ + I particles of type 1 at dual distance 
1 of each other. A pillar is a particle of type 1 at dual distance 1 of the beam not located at one of 
the two ends of the beam. The particle in the beam sitting next to the pillar divides the beam into two 
sections. The lengths of these two sections are > and sum up to £. The support of a pillared beam is 
the union of all the tile supports. The support consists of three rows ( or columns ) of sites - an upper, 
middle and lower row (or column) - which are referred to as roof, center and basement fsee Fig. 1131) . 



Figure 13: A south-pillared horizontal beam of length 10 with a west-section of length 4 and an 
east-section of length 6. 

Note that a beam can have more than one pillar. Lemma |5. II implies the following. 

Corollary 5.3 Let rj be a configuration containing a pillared beam b such that supp{b) is not 2-tiled. 
Then the configuration rj' obtained from rj by 2-tiling supp{b) satisfies H{ri') < H(q). 

5.2 Six energy reduction mechanisms 
5.2.1 Moving unit holes inside 2 tiled clusters 

In this section we show how a unit hole can move inside a 2-tiled cluster. In particular, we show that 
such motion is possible within an energy barrier QU by changing the configuration only locally. 

Definition 5.4 A set of sites S inside A obtained from a 4 x 4 square after removing the four comer 
sites is called a slot. 

Given a slot 5, we assign a label to each of the 12 sites in S as in Fig.[T3](a): first clockwise in the center 
of S and then clockwise on the boundary of S. We call the pairs 53) and (5*2, 5*4) slot-conjugate 
sites. 

Lemma 5.5 Let S be a slot, and let rjo be any configuration such that all particles in S have the same 
parity. W.l.o.g. this parity may be taken to be even, so that ri{Si) = and rj^S^) = 2. Let rji be the 
configuration obtained from rj by interchanging the states of Si and S3. Then H^rjo) — H{rii), and 
there exists a path u: rjo rji that never exceeds the energy level i?(?7o) + 6C/. 

Proof. W.l.o.g. we take rjo as in Fig. [TW b-c). Let a —?' b denote the motion of a particle from site a 
to site b. For the path uj we choose the following sequence of moves: »S'4 — y S^] S3 — y •5*4 5 S2 — ^ •5*3 5 
Si —> S2; S4 Si; S3 —y S4. The first three moves and the second three moves each are a rotation by 
■| of the subconfiguration at the sites Si, S2, S3, S4. Note that all configurations in uj have the same 
number of particles of each type and hence the changes in energy only depend on the change in the 



18 



















5*5 












Si 


S2 


S7 






Su 


S4 












Sio 


S9 



















(a) 
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Figure 14: (a) labelling of the sites in the slot (standard representation); (b) example of rjQ in the 
slot (standard representation); (c) example of rjo in the slot (dual representation), (d) rji in the slot 
(standard representation); (e) of 771 in the slot (dual representation). 



number of active bonds. Let Map be the loss of the number of active bonds between the rotating 
particles and the fixed particles, and M/j the loss of the number of active bonds between the rotating 
particles. We must show that M^p + < 6 during the six moves. To that end, we first observe that 
Mjip < 6, since the total number of active bonds between the rotating particles and the fixed particles 
is at most 6 (see Fig. HH^b)), and that M^p = 6 only after the first three moves are completed, i.e., 
when the configuration is such that all the rotating particles have a different parity with respect to the 
parity they had in configuration rjQ (recall that particles with different parity cannot share a bond). 
Next we observe that, by the choice of oj, the value of can only be or 1, and that Mn — after 
the first three moves are completed. ■ 

Lemma implies the following. 

Proposition 5.6 Let rj be a 2-tiled configuration with a unit hole. Then the configuration 77' obtained 
from 77 by moving the unit hole elsewhere satisfies H{ri') = H{ri) and (^[ri^if) < H{ri) + 6U . 

A possible 6C/-path for a unit hole inside a 2-tiled cluster is given in Fig. [131 This path is obtained 
through an iteration of local moves as explained in Fig. UM 

m 

m 
a 



-<i 



Figure 15: Motion of a unit hole inside a 2-tiled cluster. 




5.2.2 Adding and removing 12 bars from lattice-connecting rectangles 

Lemma 5.7 Let rj be a configuration consisting of a single 2-tiled lattice- connecting rectangle. Then 
the configuration rj' obtained from r/ by, respectively, 

1. adding a 12 -bar of length i > £* , 

2. adding a 12-bar of length £ < £* , 

3. removing a 12-bar of length £ > £* , 

4. removing a 12-bar of length £ < £* , 
satisfies, respectively. 
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1. H{rf) < H{rj) and <P{r], r]') < H{r]) + 2Ai + 2A2 - 4f/, 

2. H{ri') > H{r]) and <Z'(r/, rj') < H{r]) + 2Ai + 2A2 - 4U , 

3. H{r]') > H{ff} and ^{tj, r]') < H{rj) + {£ - 2)e + 4C/ - Ai, 

4. H{f]') < H{r]) and ${rj,r]') < H{rj) + {£ - 2)e + 4U - Ai. 

Proof. Recall the computations in Sections 13.11 and 14.11 

Adding a 12— bar. Adding a 12-bar of length ^ on a lattice-connecting side of a 2-tiled rectangle 
(i.e., a side such that all the particles of type 1 on that side are lattice-connecting) can be done in two 
steps: (i) initiate the 12-bar by adding a 2-tiled protuberance (see Fig. [T5)) : (ii) complete the 12-bar 
by adding a 2-tile (in a "corner") £ — 1 times (see Fig. [T7| . This can be achieved within energy barrier 
AH = 2Ai -|- 2A2 — 4J7 by following the same moves as the reference path uj* described in Section HTT] 
The energy difference due to the extra 12-bar of length i is AH{1) = Ai — le, which changes sign at 




AH = 2Ai + A2 -4(7 > 



Figure 16: A 2-tiled protuberance is added to a side of a dual rectangle within energy barrier A2. 




AH = Ai + A2 - 4{/ < 



Figure 17: A 2-tile is added in a corner between 2-tiles within a energy barrier A2. 

Removing a 12— bar. Removing a 12-bar of length £ from a lattice-connecting rectangle can be 
done by following the reverse of the path used to add a 12-bar: (i) remove £ — 1 times a 2-tile 
from a bar; (ii) remove the last 2-tiled protuberance. This can be achieved within energy barrier 
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AH{£) ~ {£ — 2)e + AU — Ai. If the cluster consists of one 12-bar only, then the path just described 
leaves ^ + 1 free particles of type 1 inside A, which can be removed (free of energy cost) afterwards. ■ 

We use Lemma 15.71 to build a northern rectangle on top of a 12-bar as follows. 

Definition 5.8 Let h denote the vertical coordinate of the sites lying on the north-side o/9^A^. For 
a given 2-tiled rectangle r in A^, let denote the vertical coordinate of the northern-most particles of 
type lin r. Then r is said to be touching the north-side of d~A~ if br — b or b^ — b —■ ^. 

In words, a 2-tiled rectangle is said to be touching the north-side of hr if it is not possible to add 
a 12-bar on the north-side within A~. Rectangles touching the south-, east- or west-side of A^ are 
defined similarly. 

Let 5 be a horizontal 12-bar of length i.e., a 2-tiled Ixl rectangle. Suppose that all sites above b 
are vacant. Then it is possible to successively add horizontal 12-bars, say m in total, on top of b until 
the north side of the rectangle grown in this way touches the north-side of A~ . The 2-tiled rectangle 
with m + 1 rows and I columns such that b is its lower- most horizontal 12-bar is denoted by □ (6) and 
is called the northern rectangle of b. 

Lemma |5 . 71 implies the following. 

Proposition 5.9 Let rj be a configuration containing a horizontal 12-bar b of length £ > i* . Then the 
configuration rj' obtained from rj by building □ [oj satisfies H{rj') < H{rj) and 'P{ri, rj') < H{rj) -\- 2Ai -|- 
2A2 - W. 

5.2.3 Changing bridges into 12 bars 

Definition 5.10 A ( south- jbridge b consists of a beam b and two (south- )pillars at the outer-most 
sites of the (south- )basement ofb. The (south-) support of b coincides with the (south-) support ofb. If 
each of the central sites of the tiles of the ( south- ) support of the bridge is occupied by a particle of type 
2, then the bridge is said to be stable (see Fig. 1181) . 

Clearly, a 12-bar is a stable bridge. North-, east- and west-bridges are defined in a similar way. 




Figure 18: A stable bridge of length 6. 

Given a bridge b. let b denote the 12-bar obtained by 2-tiling b. Lemma [511] implies the following. 

Lemma 5.11 Let rj be a configuration containing a bridge b whose support is not 2-tiled. Then the 
configuration rj' obtained from rj by changing b to b satisfies H(rj') < H(tj). 

Lemma [5.111 leads us to the following. 

Proposition 5.12 Let rj be a configuration containing a (south- )bridge b whose (south- ) support is not 
2-tiled such that the particles of its beam are lattice- connecting. Then the configuration rj' obtained 
from rj by 2-tiling supp(b) satisfies H{rj') < H(jj) and rj') < H{rj) + 4[/ + Ai. 

Proof. Let the (south-)bridge b have length £. Label the £ -\- 1 sites of its (south-)basement as 
So, si, . . . , si, from the left to the right. In order to show that supp(6) can be 2-tiled within energy 
barrier 4J7 -|- Ai , it is enough to show that within the same energy barrier a particle of type 1 can be 
brought to a site of the basement of b (from the left) that is empty or is occupied by a particle of type 
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2. W.l.o.g. si may be assumed to be such a site. The configuration thus obtained has an energy that 
is at most the energy of the original configuration (see Lemma 15. ip . The claim follows by noting that 
the particles of type 1 at the extremal sites si and Si are the two pillars of a (south-)bridge of length 
£ — I whose basement consists of the sites si, S2, ■ ■ ■ , si. 

It remains to show how a particle of type 1 can be brought to site si. Label the site north-west of 
si by vi , and the site north-east of vi by as V2- Two cases need to be distinguished: 

(1) If ?7(si) = 0, then, by the same argument as in the proof of Lemma l5.5l it is easy to show that the 
particle of type 1 at V2 can be moved to si (to obtain a configuration fj with H{fj) < H{ri)) without 
exceeding energy level H{ri) + AU. The configuration rj' is reached within an energy barrier Ai by 
bringing a particle of type 1 inside A and moving it to V2. 

(2) If r]{si) = 2, then consider the following path. First detach {AH = 2U) and remove (AH = — Ai) 
the particle of type 1 at V2, and afterwards detach {AH = 2U) and remove {AH = — A2) the particle 
of type 2 at 113. Next, move the particle of type 2 at site si to site vi {AH < 0; this particle has at 
most 2 active bonds when it sits at si), and finally bring a particle of type 1 {AH = Ai) to site V2 
{AH = —2U). Call this configuration fj. Note that H{fj) < H{ri), since effectively a particle of type 
2 with at most two active bonds has been removed, and <P{ri,if) = H{rj) + AU + Ai. Finally, observe 
that 77' is the same configuration as 77 in Case (1). ■ 

5.2.4 Maximally expanding 2 tiled rectangles 

The mechanism presented in this section, which is called north maximal expansion of a 2-tiled rectangle, 
is such that it can be applied to a 2-tiled rectangle whose north-side is lattice-connecting (even though 
this condition is not restrictive). South, east and west maximal expansion of a 2-tiled cluster are 
analogous. 

Definition 5.13 The north maximal expansion comes in two phases: a growing phase and a smoothing 
phase. 

(i) The growing phase consists of the following three steps repeated cyclically: 

1. If the particles of type 1 on the south-side of the rectangle, either at the beginning or obtained 
after step\3i constitute a south-pillared beam bs, then change supp{bs) into a 12-bar. 

2. If the particles of type 1 on the east-side of the rectangle, obtained after step [7J constitute an 
east-pillared beam h^, then change supp{he) into a \2-har. 

3. If the particles of type 1 on the west -side of the rectangle, obtained after step constitute a 
west-pillared beam b^, then change supp{bm) into a 12-bar. 

The growing phase ends after three consecutive steps leave the configuration unchanged. 

(ii) The smoothing phase consists of removing all the particles of type 2 that are adjacent to the ones 
on the sides of the rectangle that is built during the growing phase. Note that these particles have at 
most two active bonds (otherwise it would be possible to identify another pillared beam), and therefore 
removal of these particles lowers the energy. 

The outcome of the north maximal expansion (see Fig. [TO)) of a 2-tiled rectangle is again a 2-tiled 
rectangle, containing the old rectangle and such that the northern-most 12-bar of the new rectangle 
has the same vertical coordinate. 

Given a 2-tiled rectangle r, let Ti7 (r) denote the north maximal expansion of r. Corollary 15.31 
implies the following. 

Lemma 5.14 Let rj be a configuration containing a 2-tiled rectangle. Then the configuration rj' ob- 
tained from rj via (north) maximal expansion of this 2-tiled rectangle satisfies then H{ri') < H{ri). 
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(d) (e) (f) 

Figure 19: Example of north maximal expansion of a 2-tiled rectangle. The outcome of the steps 
of the growing phase are represented in pictures (b-e) , while the outcome of the smoothing phase is 
represented in picture (f). 



Lemma [5. 141 leads us to the following. 

Proposition 5.15 Let rj he a configuration containing a 2-tiled rectangle r whose north- side is lattice- 
connecting. Then the configuration rj' obtained from rj after replacing r by (r) satisfies Hij]') < H(r]) 
and <P{7], ri') < H{t]) + lOU - Ai. 

Proof. If TZ^ (r) = r, then there is nothing to prove. Therefore suppose that r is such that one its sides 
is a pillared beam. W.l.o.g. we may assume that the south-side of r is a beam b with a south-pillar. We 
must show that the south-support of b can be turned into a 12-bar within energy barrier IOC/ — Ai. 

Since supp(6) is not a 12-bar, a pillar can be chosen in such a way that at least one of the 2~tiles 
of the support the pillar belongs to (i.e., the first tile of each section of the support, counting from the 
pillar) is not a 2-tile. W.l.o.g. we let this tile be the first tile of the right-section and call it t. Let v 
denote the tile adjacent to the right site of v. In the following, the term superficial refers to tiles that 
are in the top tile-bar of the rectangle. In analogy with the proof of Lemma 15.11 several cases need to 
be considered (we stick to the order in Fig. [T^ . 

(i) {r]{tc),r]{te)) — (0,0). A particle of type 2 has to be brought to site t^ and a particle of type 1 to 
site tc- First bring a particle of type 2 to site tg, to reach a configuration i), and then proceed as 
in Case (jn]). As we will see in Case ([n]), since H{f}) = H{ri) — 3U -t- A2, the second part of the 
path can be completed without exceeding energy level H{ri) + 6U -\- A2. To reach configuration 
■q, move the particle of type 2 of the 2-tile above t to site t^ to reach a configuration called 77'. 
This can be done without exceeding energy level H{ri) 6U. Note that H{ri') = H{ri) + U . The 
unit hole that has been created at the central site of the tile above t has to be filled. This can be 
done (see Lemma 15. 5p by first moving the unit hole until it becomes superficial (configuration fj 
with energy H{fi) = H{rj')) without exceeding energy level H{ri') -t- 6C/, and then filling this unit 
hole with a particle of type 2 within energy level H{jf) -t- [/ — Ai -f A2 = H(rj) + 2U — Ai + A2. 
Thus, 77' can be reached without exceeding energy barrier 6U + A2. 
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(ii) {ri{tc),'ri{te)) = (0,2). A particle of type 1 has to be brought to site tc- Depending on the state 
of site Ve, there are three cases. 

(a) Site Ve is occupied by a particle of type 2. Move the particle of type 1 at site ti, to site 
tci to reach a configuration r/ with energy H{ri') < H{ri) + 2U within an energy barrier of 
6U. The vacancy at site ti, can be moved (again by Lemma I5.5P to the north-side of the 
rectangle within energy barrier 6C/, to reach a configuration 77 with H{fi) < H(r]), and then 
filled with an extra particle of type 1. Thus, rj' can be reached without exceeding energy 
level H{ri) + 8U. 

(b) Site Ve is empty. Move the particle of type 1 at site tj, to site Ve {AH < 3U), and then to 
site td {AH — 0). Call this configuration rj' , and note that H{rj') < H{rj) + 2U . Arguing 
as above, we see that the vacancy at site tb can be filled without exceeding the energy level 
H{-q)+m. 

(c) Site Ve is occupied by a particle of type 1. Observe that the particle of type 1 at tb has fc < 3 
active bonds and the particle of type 2 at Ve has m < 2 active bonds. It is possible to move 
the particle at site Ve to site tc {AH — (m — k)U), and then the particle at site to site 
Vc {AH — {k — m)U). The configuration 77', reached within energy barrier (fc — m)U , has 
energy H{r]') < H{rj) + kU . Again, the vacancy at site if, has to be filled with a particle of 
type 1. This can be done without exceeding the energy level H{rf) + (6 + k)U < H{rf) + QU . 



(iii) {ri{tc),'i]{te)) = (2,0). The particle of type 2 at site tc is moved to site te without increasing the 
energy. Then argue as in Case (ju]). 

(iv) {r]{tc),r]{te)) — (2,2). The particle of type 2 at site tc has to be replaced by a particle of type 
1. Remove the particle of type 2 aX te- To do this, first create a superficial unit hole (which can 
be done within energy barrier AU — Ai by creating a hole in a corner tile of the rectangle) and 
move this vacancy to site te- By Lemma l5.5l this can be achieved without exceeding energy level 
H{riQ) + 10[/ — A2. Then move the particle of type 2 at site tc to site te {AH < 0). Call 77' the 
configuration that is reached in this way. Note that H{rj') < H{rj) — A2 + 'iU . To bring a particle 
of type 1 to site t^ argue as in Case (ju]), to arrive at H{ff) < H{rf) + 12U — A2. 

(v) {ri{tc),'ri{te)) = (1,0). A particle of type 2 has to be brought to site te- Move the unit hole at te 
to the top tile-bar of the rectangle. This does not change the energy of the configuration and can 
be done within energy barrier 6U by Proposition 15. 61 The task reduces to filling a superficial unit 
hole on the surface of the cluster with a particle of type 2. This can be achieved within energy 
barrier U + A2 — Ai. Therefore the maximal energy level reached in this case is H{ri) + 6U. 

(vi) {ri{tc),'ri{te)) — (0,1). Move the particle of type 2 from site te to site tc- This move does not 
increase the energy of the configuration. Then proceed as in Case (jvj) . 

(vii) {ri{tc) , ri{te)) = (2,1). The occupation numbers of sites tc and te have to be exchanged. To 
do this, first remove the particle of type 1 at site tb to obtain a configuration rj' with energy 
H{rj') < H{rj) + 3[/ without exceeding the energy level H{rf) + 10[/ — Ai (again use Lemma [^31) . 
Move the particle of type 1 from te to tb {AH < 0) and the particle of type 2 from tc to te 
{AH = 0). Call i) the configuration that is reached in this way. Note that H{ff) < H{rj) + [/ — Ai. 
Proceed as in Case ([u]) to conclude within energy barrier of IOC/ — Ai. 

(viii) {ri{tc), vi^e)) = (1, 1)- The particle of type 1 at site te has to be replaced by a particle of type 2. 
This can be done as follows. First the particle of type 1 sitting a site tb is removed. To achieve 
this, first remove a particle of type 1 at the north-side of the rectangle and then (use Lemma [575)1 
move the vacancy to site tb- The configuration that is reached, which we call 77', is such that 
H{ri') < H{r]) + 3/7 — Ai. Next, move the particle of type 1 at te to site tb {AH — 0), to reach 
a configuration 77 whose energy is H{fi) = H{rf) — Ai. Finally, argue as in Case (jvj), to arrive at 
H{fi) <H{rj)+3U-Ai- 
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Finally, note that PTT^ implies max{6C/ + As, IOC/ - Ai, 12U - A2} = IOC/ - Ai. 

By Lemma [5. 11 H{r}') < H{'ri), and therefore the same argument can be used to show that all the 
right-sections of the support can be 2-tiled within the same energy barrier. The left-section can be 
2-tiled analogously. 

To conclude, it remains to be shown how particles of type 2, possibly adjacent to one side of the 
rectangle, can be removed from A. Call t the tile associated with the particle p of type 2 that has 
to be removed {p sits at site te) and v the tile adjacent to t belonging to the rectangle. First bring 
a vacancy to site Ve within energy barrier IOC/ — A2 (one way to achieve this has been described in 
Case ((iv|) above) and then move p to site Ve (see Lemma [575)) . ■ 



5.2.5 Merging adjacent 2 tiled rectangles 

Definition 5.16 A 12-bar bi of length £ of a cluster ci is said to be adjacent to a 12-bar 62 of length 
m < £ of a cluster C2 if there exist m mutually disjoint pairs (gj, q\) of particles of type 1 with q\ S 61 
and 52 € &2 such that u(c[\) — ^(92) = with \\v\\ = \\f2 for i = 1, . . . , m. The vector v is called the 
offset of 62 with respect to bi . The tiles in bi have a different parity than the tiles in &2 • The particles 
q{ G bi, i — 1, . . .m, are called the external particles of bi with respect to 62; o,nd the particles q\ G &2> 
i — 1, . . . , m, are called the external particles of &2 with respect to bi . 

Proposition 5.17 Let rj be a configuration that contains two adjacent 2-tiled rectangles. Then the 
configuration rj' obtained by "merging" these two rectangles satisfies H[ri') = H{ri) and <P{rj, rj') < 
H{r,) + 2U-Ai. 

Proof. Given two adjacent bars bi and 62 with offset v = {vi,V2) in a configuration 77, we want to 
define the sliding of 62 onto bi along v. The resulting configuration 77' is such that all the particles 
of type 2 originally in 62 are slid by (ui,W2) with respect to their position in rj, and all the external 
particles of type 1 of 62 with respect to bi are slid by (wi, —V2) when the two bars are horizontal and 
hy(—vi,V2) when the two bars are vertical. Via the sliding, the m 2-tiles in 62 are turned into m 
2-tiles with the same parity as the tiles in bi. It is easy to see that H{ri') = H{r]), since neither the 
total number of active bonds of the configuration nor the number of particles of each type is changed. 

To describe the sliding of a bar onto another bar along a vector v, we may assume w.l.o.g. that 
the two bars are vertical and that the vector v is equal to (— |, — ^) (Fig. 20(a) I. Start by moving the 



lower-most external particle of type 1 in 62 over the vector v' = (|, — i) (Fig. 20(b) I. This leads to an 
increase by U in energy. Then move the lower- most particle of type 2 over the vector v (Fig. 20(c) ). 
Since the number of deactivated bonds is equal to the number of new bonds activated, this move does 
not change the energy. Proceed by moving over the vector v' the second particle of type 1 from the 
bottom of the bar (Fig. |20(d)[ ). This also is a move that does not change the energy. Afterwards, the 



second particle of type 2 from the top is moved over the vector v (Fig. 20(e) I. This sequence of moves 
proceeds iteratively (without a change in energy) until the TO-th particle of type 2 has been moved over 
the vector v. Finally, the {m + l)-st external particle of type 1 is moved over the vector v' (Fig. [20(f)] ). 
This move decreases the energy by U. Thus, U is the energy barrier that must be overcome in order 
to realize the sliding of a 12-bar onto another 12-bar over the vector v. 

It is clear that, given a configuration 77 containing two 2-tiled rectangles ci (with vertical side length 
£) and C2 (with vertical side length m < £) with offset v, it is possible to reduce 77 to a configuration 
rj' such that ci and C2 are merged into of a single cluster by sliding one bar after another, without 
exceeding energy barrier AiJ = U , provided the other clusters of rj do not interfere with this procedure. 
Sliding the last bar of C2 we get an excess of free particles of type 1, which can be removed from A, 
lowering the energy. In particular, the configuration rj' obtained via the sliding of C2 onto ci along 
v without exceeding energy level Hij]) -\- U has energy H{r]') = H{ti) — {m -\- l)Ai, since the two 
configurations consist of the same number of 2-tiles, and 77' contains 771 -|- 1 particles of type 1 less than 
77. Moreover, <P{rj,if) = H{rj) + U . 

In the argument above, the first move consisted of moving down-right a particle of type 1 of 62 
to an empty site (say, site i). If in configuration 77 site i is occupied by a particle of type 1, then 
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(d) (e) (f) r,i 

Figure 20: The sliding of 62 onto bi. 



the shding of the vertical 12-bar can be realized by modifying the procedure as follows. First remove 
from the box the top-left particle of type 1 of 62 sitting at site j to reach a configuration with energy 
H{ri) + U — Ai (which can be done without exceeding energy level H{ri) + U). Then move to j the 
particle of type 1 sitting at site /c = j + d = j + (— ^, — ^) in 77, which increases the energy up to level 
Hijf) + 2U — Ai. Then site k is filled with the particle of type 1 originally at site k + (i, — i) without 
an increase in energy. It is possible to continue in this way until the configuration obtained after the 
first step of the above case is reached. This configuration has energy H{rj) + J7 — Ai. Then proceed as 
in the above case until 62 is slid onto 61. This leads to a configuration with energy H{rj) — Ai < H{rf). 
In order to perform the (modified) sliding procedure, it is sufficient to assume that the north-side of 
rectangle C2 is lattice-connecting. ■ 



5.2.6 Removing subcritical clusters 

The cleaning mechanism defined in this section produces a configuration for which we have a certain 
control on the geometry of the constituent clusters. In particular, these clusters will be suitable for the 
application of the previous five energy reduction mechanisms. We begin by looking at pending dimers 
(see Fig.HI]). 
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Figure 21: A pending dimer is the pair of particles circled in the picture. 



Definition 5.18 A pending dimer consists of two adjacent particles of different type such that the 
particle of type I is lattice- connecting and has only one active bond and the particle of type 2 has at 
most three active bonds. 

Proposition 5.19 Let rj be a configuration containing pending dimers. Then there exists a configura- 
tion T]' not containing pending dimers that satisfies H{if) < H{ri) and <P{r],ri') < H{r]) -f- SU -f A2. 
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Proof. If the particle of type 2 has at most two active bonds, then simply remove the pending dimer. 
This reduces the energy, since two bonds are deactivated and a particle of each type is removed from A 
(AH < 2U — Ai — A2 < 0), and can be achieved within an energy barrier 2U — Ai along the following 
path: first detach [AH — U) and remove {AH — — Ai) the particle of type 1, then detach {AH < U) 
and remove {AH = — A2) the particle of type 2. 

If the particle of type 2 has three active bonds we have two cases: 

(i) The fourth neighbor of the particle of type 2 of the pending dimer is empty. In this case rj' is 
obtained by filling this empty site with a particle of type 1 in order to obtain a 2-tile, which 
lowers the energy since Ai < U. To do this, temporarily remove the pending dimer as described 
above. This leads to a configuration ij with energy H{fj) = H{rj) + 3C/— Ai — A2 reached within 
energy barrier 3C/ — Ai. Then bring a particle of type 1 to the designated site {AH < Ai) and 
finally put back the dimer. The whole path is realized within energy barrier 3U + A2. 

(ii) The fourth neighbor of the particle of type 2 is occupied by a particle of type 2. In this case 
rj' is the configuration such that the dimer is removed and the site originally occupied by the 
particle of type 2 of the dimer is occupied by a particle of type 1. To obtain 77' from 77, remove 
the pending dimer (again, as above, within energy barrier 3U — Ai), to reach a configuration 77 
with energy H{Ti = H{'q) + 3[/ — Ai — A2, and bring a particle of type 1 within energy barrier 
Ai. To conclude, observe that ^(77') = H{ri) + 2?/ - A2 < H{ini). 



The cleaning mechanism works as follows: 

1. Remove all the lattice-connecting free particles from the configuration. 
After that repeat cyclically the following two steps: 

2. Iteratively remove/transform all the lattice-connecting pending dimers. 

3. Bring a particle of type 1 to any of the free sites adjacent to the lattice-connecting particles of 
type 2. 

Repeat the cleaning mechanism until the configuration is not affected anymore. Each of the three steps 
can be performed within energy barrier 3C/ -|- A2. Moreover, each step reduces the energy. 

Lemma 5.20 The outcome of the cleaning mechanism is either a configuration such that the first 
particle encountered while scanning A in the lexicographic order is a particle of type 1 belonging to a 
horizontal stable (south-)bridge, or the configuration □. 

Proof. Call q the first particle of A in the lexicographic order. Recall that the dual coordinates of q 
are denoted by u{q) — 1*2(9)). Step[3]of the cleaning mechanism guarantees that g is a particle 

of type 1. The fact that q is the first particle in the lexicographic order implies that: (i) all the sites 
above u{q) are empty; (ii) all the sites with the same vertical coordinate as q lying on the left of q are 
empty as well. As a consequence of (ii) , all the sites on the left of q with vertical coordinate {q) — \ 
are lattice-connecting and therefore cannot be occupied by a particle of type 2. Since q cannot be a 
free particle, the site with coordinates {ui{q) + i, U2{q) — \) must be occupied by a particle p of type 
2. Let s{p) be the longest sequence of tiles adjacent to t{p) such that the central site is occupied by 
a particle of type 2. Obviously, p is the left-most particle of type 2 in s{p). Call p the last particle of 
type 2 in s{p) and q the particle of type 1 with coordinates {u\{p) -\- ^, M2(p) -I- \). (Note that p and p 
may coincide.) All the sites on the north-side of s{p) are lattice-connecting and hence are occupied by 
a particle of type 1. To conclude, observe that both p and p must be saturated, otherwise at least one 
of the pairs (g,p) and {q,p) constitutes a pending dimer. ■ 
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5.3 Energy reduction of a general configuration: Proof of Theorem 11.41 



Fix any rj ^ {□, ffl}. In this section we will give a general procedure, called energy reduction algorithm, 
that allows us to construct a path a;: r/ — !• r/r with rjr G {□, ffl} such that niax^g^ H{£^) < H{r]) + V* 
with V* < lOU - Ai and H{r]r) < H{ri). Note that if r/^ = ffl, then H{r]r) < H{r]) because A^stab — ™- 
The construction uses the six energy reduction mechanisms described in Sections I5.2.1H5TT^ and relies 
on Propositions 15.61 15.91 15.121 15.151 15.171 15.191 which are the key results of these sections. The 
maximal energy barrier in these propositions is IOC/ — Ai. Note: The energy reduction mechanisms 
in Sections 15.2.21 and 15.2.31 concern single droplets far away from d~A and have an energy barrier not 
exceeding 4C7 + Ai < F* (see below (|1.14p ). For such configurations, the energy can be essentially 
reduced by saturating particles of type 2 and by adding and removing 12-bars. This explains the 
remark made in Section [L^ item 4. 

In the remainder of this section we call supercritical a 12-bar of length > i* . Similarly, we call 
supercritical a dual rectangle with both side lengths > £* . 

Proof. As a preliminary step, perform the cleaning mechanism. If the outcome is □, then the claim 
is proven. Otherwise, let bi be the first bridge encountered in the lexicographic order (which exists by 
Lemma r5.20p . This bridge can be turned into an 12-bar bi (see Section [5.2.31) . If the length of bi is 

< £*, then the 12-bar bi can be removed, which lowers the energy (see Section [5221) • this case, go 
back to performing the cleaning mechanism. W.l.o.g. we may therefore assume that the length of bi is 
> £*. 

By construction, all sites above bi are empty, and therefore it is possible first to construct the 
2-tiled rectangle ri — n within energy barrier 2Ai + 2A2 — 4L/ (again lowering the energy), and 
then expand ri to the rectangle Ri — TZ^ (ri) (see Section [5.2.41) . If the vertical side length of i?i is 

< i* , then Ri can be removed (lowering the energy), and it is possible to perform again the cleaning 
mechanism. 

Therefore suppose that Ri has both its side lengths > i* . In the remainder of the section we 
will show how to reach within energy barrier lOU — Ai a configuration containing a rectangle Rnw 
touching both the north-side and the west-side of A~ whose support contains the support of Ri . Once 
this has been achieved, it is possible to argue for Rnw in the same way as for Ri in order to reach 
a configuration containing a rectangle Rnwe touching the north-side, the east-side and the west-side 
of A~ whose support contains the support of Rnw- Repeating the same argument for Rnwe, it is 
possible to reach ffl. 

The construction of Rnw is obtained by using an algorithm called invasion of Ri, which is con- 
structed with the help of techniques similar to the ones that were used to build 

(A) Invasion of See Fig. [22] Let (ai,6i) be, respectively, the horizontal and the vertical 

coordinate of the left lower-most particle of Ri (which is of type 1). Define A(i?i) C A to be the set 
consisting of the sites whose vertical coordinate is > 61 and horizontal coordinate is < ai. In words, 
A(i?i) contains the sites of A on the left of Ri. Perform the cleaning mechanism (see Section [5.2.6^ 
and scan A(i?i) in the lexicographic order. Three cases are possible. 

1. A(i?i) is empty. Add, if possible (i?i might already be touching the west-boundary of A~), 12-bars 
onto the left side of Ri until the resulting cluster touches the west-boundary of A~ . 

2. The first horizontal bridge 62 encountered in A(i?i) has length < £*. Remove the particles of the 
(south)-support of the bridge, lowering the energy of the configuration, and restart the covering of 
A{Ri). 

3. The first horizontal bridge 62 encountered in A(i?i) has length > £*. As for 61, first turn 62 into the 
12-bar &2, then build the 2-tiled rectangle r2 = n (62), after that expand r2 to i?2 — TV (r2), and 
finally perform the cleaning mechanism. Note that the support of R2 may cover (part or possibly all 
of) the support of This means that during the maximal expansion, some of the sites of supp(i?i) 
were in the support of the pillared beam that is going to be 2-tiled. Each time this happens, R2 
absorbs an entire vertical supercritical 12-bar of i?i (see Section [5. 2. 4p . Call i?i what is left of i?i 
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after the maximal expansion of i?2 • The following three cases are possible: (i) Ri does not contain 
any particle (i?i = 0); (ii) Ri -< Ri (in the proper sense); (iii) Ri = Ri. In Case (ii), the rectangles 
i?2 and i?i are necessarily adjacent (more precisely, the right-most 12-bar of i?2 is adjacent to the 
left-most 12-bar of whereas in Case (iii) the two rectangles may or may not be adjacent. Note 
that this implies that if Ri ^ Ri, then R2 is necessarily supercritical. Obviously, if i?i 7^ 0, then it 
is again a 2-tiled rectangle, and there are several possibilities. 

(a) i?2 is not supercritical. This implies that i?i — Ri. Remove R2 from A, put i?i = Ri and 
restart the invasion of i?i . 

(b) i?2 is supercritical and i?i = 0. Change the name of i?2 to i?i and restart the covering of A(i?i). 

(c) i?2 is supercritical and is adjacent to Note that both rectangles touch the north-side of 
A^. Call 7?™^=^ the rectangle with the largest vertical length (in case of a tie, w.l.o.g. choose 
i?i) and call i?™'" the other rectangle. Slide i?™'" onto i?™'*''. This is possible because the 
smoothing phase of the maximal expansion (see Section I5.2.4P removes all the particles of type 
2 that may interfere with the sliding of the 12-bars. Then perform again the maximal expansion 
of R"^^^, i.e., the rectangle that has not been moved during the sliding. These steps bring the 
configuration to a rectangle whose support contains supp(i?2) U supp(i?i) U A(i?i). Call this 
rectangle i?i and restart the invasion of 

(d) i?2 is supercritical and is not adjacent to Ri. This implies Ri = Ri. Start the invasion of R2 
(see below). 

In order to complete the proof, it remains to show how the invasion of R2 carries over. To that 
end, we introduce the following recursive algorithm realizing the invasion of Ri for i = 2,3, . . ., etc. 

(B) Invasion of Ri. Call Ri-i what is left of after the invasion of Ri+i. There are three cases: 

I. Ri-i = (i.e., the support of Ri-i is completely covered by Ri). Put Ri^i ~ Ri and restart the 
invasion of Ri-i- 

II. Ri-i 7^ and Ri and Ri-i are adjacent. Call R'°^^^ the rectangle with the largest vertical side 
between Ri and Ri-i (in case of a tie, w.l.o.g. choose i?™^'^ = R^) and call i?™'" the other 
rectangle. Slide i?""" onto R'^^^^ and perform the maximal expansion of R"^^^. Call Ri-i the 
outcome of the maximal expansion of R^^^^ and restart the invasion of Ri-i- 

III. Ri-i and Ri and Ri-i are not adjacent. If Ri is on the left of Ri-i, then let (oj, bi) denote, 
respectively, the horizontal and the vertical coordinate of the lower right-most particle (which is 
of type 1) of Ri, and call A{Ri) the subset of A(_Rj_i) consisting of those sites whose vertical 
coordinates are > bi and whose horizontal coordinates are > a^. If Ri is on the right of then 
let (ai^bi) denote, respectively, the horizontal and the vertical coordinate of the lower left- most 
particle (which is of type 1) of Ri, and call A(i?i) the subset of A(_Ri_i) consisting of those sites 
whose vertical coordinates are > bi and whose horizontal coordinates are < a^. In words, A(i?i) 
consists of those sites of A(i?i_i) between Ri^i and Ri. Perform the cleaning mechanism and 
scan A(i?j) in the lexicographic order. There are again several cases. 

1. A{Ri) is empty. Call the rectangle with the largest vertical side between Ri and Ri-i (in 
case of tie, w.l.o.g. choose R''^^^ = i?j) and call i?™'" the other rectangle. Add vertical 12-bars 
on the side of i?™'" facing R''^^^ until (depending on the parity of the rectangles) it becomes 
adjacent (different parity) to R™'^^ or it is at distance 1 (same parity) from i?™^''. In the first 
case, slide the extended i?™'" onto R^^^. Perform the maximal expansion of R™^^, and call 
Ri-i the rectangle obtained in this way, whose support contains supp(i?i) U Ri-i U A(i?i_i). 
Restart the invasion of Ri-i. 

2. The first horizontal bridge encountered in A{Ri) has length < £*. Remove the particles 
of the (south)-support of the bridge, lowering the energy of the configuration, and restart the 
invasion of Ri. 
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3. The first horizontal bridge bi+i encountered in A{B.i) has length > £*. First turn bi+i into 
the 12-bar bi+i, then build the 2-tiled rectangle r^+i = □ after that expand r.i to 

Ri+i = Ti7 (ri+i), and finally perforin the cleaning mechanism. Call Ri what is left of Rt after 
the maximal expansion of Ri+i. The following cases are possible. 

(a) is not supercritical. This implies Ri = Ri. Remove from A, put Ri = Ri, and 
restart the invasion of Ri. 

(b) is supercritical and Ri = 0. Change the name of Ri+i to Ri, and restart the invasion 
of 

(c) Ri+i is supercritical and is adjacent to Ri. Note that both rectangles touch the north-side 
of A^. Slide the rectangle with the shorter vertical length onto the other rectangle and 
perform again the maximal expansion of the rectangle that has not been moved during 
the sliding. These steps bring the configuration to a rectangle whose support contains 
supp(i?i+i) U supp(i?i) U A(Ri). Call this rectangle Ri and restart the invasion of Ri. 

(d) Ri+i is supercritical and is not adjacent to Ri. This implies Ri = Ri- Start the invasion 
of Ri+i. 

The finiteness of A ensures that the algorithm eventually terminates. ■ 
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Figure 22: Example of invasion of the dual rectangle Only the support of the relevant clusters 
are drawn and the parity of different clusters is not indicated. The set A(i?i) contains a supercritical 
bridge belonging to cluster A (Fig. 22(a) I. Growing this bridge via the construction of its northern 



rectangle and its subsequent maximal expansion leads to the supercritical rectangle R2 (Fig. 22(b) I 



Next, the invasion of A{R2) has to be performed in order to complete the invasion of The set A(i?2) 
contains a supercritical bridge belonging to cluster B, which is grown into the supercritical rectangle 



i?3 (Fig. 22(c) ). Note that R3 partly covers the support of Ri and that R3 and Ri are adjacent. The 
invasion of R2 proceeds via the invasion of R3. Since A(i?3) is empty, the invasion of R^ is carried out 
by adding 12-bars to the left-side of R3 until R2 is at dual distance 1. After that a maximal expansion 
produces a dual rectangle that covers the support of R2 (Fig. |22(d)| . The new dual rectangle R2 is 
adjacent to Ri. The two rectangles are merged and a maximal expansion gives a new rectangle Ri 
(Fig 22(e) I. Now A(_Ri) is empty and can be filled by adding 12-bars to the left-side of Ri until the 



rectangle Rnw is obtained (Fig. 22(f)) 
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